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I. INTRODUCTION 



We use an operational approach to discuss ways to measure the higher-order cross-correlations 
between optical and matter-wave fields. We pay particular attention to the fact that atomic fields 
actually consist of composite particles that can easily be separated into their basic constituents by a 
detection process such as photoionization. In the case of bosonic fields, that we specifically consider 
- - - here, this leads to the appearance in the detection signal of exchange contributions due to both 

the composite bosonic field and its individual fermionic constituents. We also show how time-gated 
' counting schemes allow to isolate specific contributions to the signal, in particular involving different 

, , orderings of the Schrodinger and Maxwell fields. 

CD 

Q 

. 

The experimental realization of Bose-Einstein condensation in dilute atomic systems has opened up exciting 
new directions of research in atomic, molecular and optical physics. Recent developments of special interest in the 
^ : context of the present paper include the first experimental demonstration of the four- wave mixing of de Broglie waves 
in a Sodium condensate [|| , the generation of dark solitons by optically imposing a phase shift on a condensate wave 
function 0,^, the realization of several types of atom lasers P-p^, and, most importantly perhaps, the demonstration 
of superradiant Rayleigh scattering in a condensate jl^ . This latter experiment is particularly significant in that it 
I is the first example of a four-wave mixing process involving dynamically two optical waves and two atomic matter 
^\ I waves. 

On I In recent work [p^, we initiated a study of the mutual coherence between optical and matter waves, discussing 
(-H . possible measurement schemes, and illustrating our results in a simple example. The goal of the present paper is to 
extend these results to the study of higher-order coherence functions. This study is motivated by several reasons: 
^ , First, it is well established in optics that higher-order detection schemes are required to test certain aspects of the 
' nonlocal predictions of quantum mechanics. p5[ | It would be useful to carry out similar experiments with massive 
I particles. Second, the generation of a quantum entanglement between optical and matter waves is of much interest in 
applications such as the optical manipulation and control of Schrodinger fields. In particular, a number of sophisticated 
techniques have been developed to generate optical waves of prescribed statistical properties. It would be desirable 
to develop methods to transfer them to matter waves, while possibly amplifying them in the process. This might 
for example open up the way to the generation of atomic Fock states of known atomic number, which are of great 
potential interest in matter-wave interferometry. ||l6| More generally, the interplay between Maxwell and Schrodinger 
waves is the cornerstone of most potential applications of atom optics, from situations such as lenses, gratings and 
mirrors where light plays a passive role to situations of dynamical coupling, including matter-wave amplifiers p^ ] 
and ultra-sensitive atom-optical detectors. These developments require one to characterize not just the statistical 
properties of the optical and matter fields individually, but also their cross-coherence properties. 

The general detection scheme that we consider consists of one or more standard photodetectors for the optical 
field, together with matter-wave detectors which operate by annihilating an atom via photoionization. As a result, 
electrons and ions are produced, and the resulting electron current is measured. The atomic detection, while on the 
surface similar to photodection, is therefore fundamentally different in that the atoms, which are composite particles, 
are not so much destroyed by the detection process as transformed into a pair of two new particles. For example, if 
the atoms under consideration are composite bosons before detection, the resulting ions are fermions. The quantum 
statistics of the resulting particles is important in those situations involving exchange terms, as we discuss in detail 
later on. 

The paper is organized as follows: Section II reviews the model of joint atomic and optical detection of Ref. p^ ] 
and establishes the notation. Section III, which is the central part of the paper, discusses the kinds of correlation 
functions that can be detected. It accounts in detail for the quantum statistics of the ions and electrons resulting from 
photoionization, including the effects of particle exchange, and illustrates how different orderings can be achieved by a 
proper time-gating of the detectors. Section IV illustrates these results in a simple atom optics example, and Section 
V is a summary and conclusion. 



1 



II. DETECTION SCHEME 



We consider the joint detection of the coherence properties of dynamicaUy coupled optical and matter-wave fields, 
the latter one being assumed to be bosonic for concreteness. In Ref. Q, we proposed a detection scheme that can 
achieve this goal. It consists of a series of standard photodetectors, which we call in the following Maxwell detectors 
and of matter-wave Schrddinger detectors. The Schrodinger detectors operate by tightly focussing a laser beam on 
the atoms, ionizing them, and measuring the resulting electron current. We describe the optical field by its electric 
field operator E(r,t) and the matter-wave field by the (multicomponent) Schrodinger field operator '4'(r, i), which 
satisfies the bosonic commutation relation 

[vl>4r),*t(r')]^^„^5(r_r'). (1) 

Here r is the center-of-mass coordinate and the indices a and (3 label the internal (electronic) state of the atoms. 
The dynamics of the compound system is governed by the Hamiltonian 

Ho = n,n + ns + v (2) 

where Tim and TLg describe the evolution of the light field and the matter-wave field, respectively, while V is the 
interaction between these two fields, typically the electric dipole interaction. 

In Ref. [Q, we considered the measurement of the lowest-order joint-correlations between the two fields. In that 
case, just one photodetector and one atom detector were required. In contrast, the measurement of higher-order 
correlations involves the use of several light and matter-wave detectors. 

As in Glauber's photodetection theory [|l9|, the photodetectors are modeled as single two- level atoms whose excited 
state is in the continuum, corresponding to photoionization. The interaction of the system with N such detectors is 
described as usual by the electric dipole interaction Hamiltonian 

N 

V„, = ^d,.E(r,) (3) 

i=i 

where dj is the electric dipole moment between the ground electronic state and the continuum state of the j^^ 
detector, located at the position r^. We assume for simplicity that the dipole moments are real and parallel to 
the field. Similarly, the interaction between the matter-wave field and the M Schrodinger detectors, assumed to be 
point-like and at locations Vj , is given by 

N+M 

Here, ^ and are the field operators associated with atoms in the (ground state) sample Q and the atomic continuum 
state respectively, 51^ is the coupling constant proportional to the Rabi frequency of the ionizing lasers of frequency 
ujh ■ As pointed out in Ref. p8| , the assumption of point-like detectors is justified for atomic samples at temperatures 
well below the recoil limit, in which case the ionizing lasers can easily be focused to a spot size much smaller than 
the thermal de Broglie wavelength. 

The atoms in the Schrodinger field '4'(r) are actually composite bosons which are photodissociated by the laser 
into an electron and an ion. Hence, the creation operator *J,(r) should actually be understood as the product of two 
fermionic creation operators. Since what is detected is the electron current, it is important to treat this property of 
the system properly. We proceed by expanding the continuum field ^ (r) in terms of plane waves of momentum q as 

*,(r) = y"rf3#q(r)6,q (5) 

where 



^We consider only scalar fields in the following for notational simplicity. The extension to the detection of multicomponent 
Schrodinger fields is straightforward. 
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Mr) - (^e-^ 



(6) 



and 



the corresponding Hamiltonian being 



uj^q^hq^/2M + Lu^, (7) 
[b..i,bl^,]=Siq-ci')S,,,, (8) 

ni = nY^ I (fiqto^qhl^K^. (9) 

We now decompose the atomic creation operators 6^^ in terms of products of electron and ion creation operators. We 
assume for concreteness that the atoms under consideration are spin-zero particles, so that the spins of the resulting 
electrons and ions are opposite. As a result of the large difference between the ion and the electron masses, the kinetic 
energy of the system is carried almost entirely by the ion. If furthermore the photon energy is insufficient to excite 
the ion above its ground electronic state, we have that 

bt^^Y.<s f d'k^t{k)cl^,{k) (10) 

where fuik) are momentum-space electron wave functions. The sum appearing in this equation is over the spin 
states s = ±1/2, and the ion and electron annihilation operators, Oq^^ and Ci,^s(k) respectively, satisfy the Fermi 
anticommutation relations 



and 



[aq,s,al,^^,]+ ^ 6syS{q~q') (11) 
[c,,,(k),4,.,,(k')]+ = ^.,.'V.'(5(k-k'). (12) 



Combining Eqs. (^ and (^) gives 

^T.J d^qMr)a^,s J d^kip,{-k)c,,-s{k). (13) 

The Schrodinger detectors measure the electron current, but are insensitive to the final state of the ions, which will 
therefore be summed over. 

The total Hamiltonian describing the coupling of the system to the detectors is finally 

Vd = V™ + V, (14) 

and the full Hamiltonian including system, detectors and interaction is 

n = na + ni + nA + Vd (is) 

where TCa is the Hamiltonian of the Maxwell detectors, modeled as usual as two level-atoms. 

III. COUNTING SIGNAL 

The probability amplitude of exciting all detectors is given to lowest order by (iV-|- A/)-th order perturbation theory. 
To that order, the probability amplitude for the transition from the initial state \i) to a final state [/) is given by 



^ , N+M i-t 



d-Tl / dT2... / dTN+M{f\Vd{Tl)Vd{T2)...Vd{TN+M)\i) (16) 

Jo Jo 
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where ti > T2 > ... > tn^i\j, and the interaction with the detectors is now expressed in the interaction picture as 



where Hd = "Hi + Ha is the free Hamiltonian of the detectors. 

When substituted into Eq. (]T6|), the interaction Hamiltonian ( p^ yields (A^ + M)*^^+^^ terms. It is well known that 
for photodetectors operating by absorption, only the positive frequency part of the electric field operator contributes 
to the signal. A similar situation also occurs for the Schrodinger detectors provided that no ion or electron is initially 
present, i.e., the field ^p{r) is initially in the vacuum state. If that is the case, the only terms contributing to 
the detection signal are those proportional to \I'J,(rj)SI'(rj)e~''^''*, which result in the detection of normally ordered 
correlation functions of the Schrodinger field. Further neglecting terms involving multiple ionization at a single 
detector reduces the probability amplitude (|l^) to 

N+M 

■ dn J dT2...J dTM+M (/|V^^(^i)vi+^(r2)...V^+_)^(rjv+M)N). (18) 

^ ^ all permutations 

In this expression, Vi can correspond to either the detection of a photon, or to the ionization of an atom. All permu- 
tations of such detection events such that a total of N photodetections of the optical field and M photoionizations 
of the Schrodinger field take place must be included in the sum. The superscrit refers to the familiar "positive 
frequency part" in the Maxwell detectors and to the ordering of the Schrodinger detectors discussed before Eq. jl^). 
Hence this expression includes {N + terms. Since the optical and atomic field operators do not generally com- 
mute at different times, it is not possible to regroup them into just one contribution. However, it is possible to use a 
time-gated detection scheme to separately detect the various correlation functions of Eq. ([isl), as we now discuss. 

We consider the experimental situation where the i-th detector is turned on at t = and turned off at some later 
time ti, so that Vd becomes 

N+M 

Mi) = [®(*) ~ - ^»)] ■ (19) 

1=1 

where Q{t) is the Heaviside function. A given ordering of the switch-off times, for example ti < t2-.. < In+m < t 
results in just one term in Eq. (|l8|) contributing to the counting rate. 

As is the case in standard optical photodetection theory , we can safely assume that the photocurrents generated 
in the Maxwell detectors are independent and distinguishable. However, more care must be taken when dealing with 
the Schrodinger detectors: Because the de Broglie wavelength of the ultracold atoms can be very large, it is not 
correct in general to assume that the electrons generated by the photoionization process are distinguishable. Hence, 
it is necessary to account for the various exchange terms arising in the detection signal jl^ . 

We illustrate how this work in the case of a detection scheme involving two Maxwell detectors at positions ri and 
r2 and turned off at times ti and t2, respectively, as well as two Schrodinger detectors at and and closed at 
times ^3 and For the time ordering ti < t2 < < < t, the only term resulting from the transition probability 
amplitude (|l8|) that contributes to the counting rate turns out to be 



^W^i E f^:,(r3)l^:,(r4)J^.3(r4)f^.Jr3) 



X 



/ du / dr^... / dTi / dT[e-''^''^^^-^^+^'-^^^ 
Jo Jo Jo Jo 

X (di(T{)di(Ti)d2(T:^)d2(r2))(*^i(r3,T^)*^,(r4,T^)*I,3(r4,r4)*I,4(i'3,T-3)) 

X (E(-)(ri,T0E(-)(r2,T^)*t(r3,r^)*t(r4,^^)*(r4,r4)*(r3,r3)E(+)(r2,T2))E(+)(ri,ri)) (20) 

where we have introduced as usual a sum over a complete set of final states. Note that in deriving this expression, 
we have factorized the expectation values into products of detector operators and of system operators correlation 
functions. This is appropriate in the absence of quantum entanglement between the states of the detectors and of the 
system to be characterized. 

In order to properly account for the quantum statistics of the detector correlation function 

Cf EE (*,,(r3,r^)*.,(r4,r^)*|;jr4,T4)*|;^(r3,T3)), (21) 
we reexpress it using Eq. ( |T3| ) as 
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X 0qi (r3)0q2 (i-4)0;3 (r4)(/>;4 (r3)¥'z.i i'ki)ip^^ 0^2)'pU 0^3)^1, (^4) 

(k4)) (22) 

where the sum is over a set of discrete and continuous indices s, q and k and we have assumed that the electron and 
ion dynamics are governed by their free Hamiltonians. Since the ion and electron fields are in the vacuum state, the 
expectation value of the ions and electrons operators can be easily evaluated using their fermion anticommutation 
relations to give 

(aqi.siOqa^s^aL.ss^L.Si) = '5(qi - q4)'5(q2 - q3)SsusJs2.s3 ~ - <\3)K^i - q4)^si,s3'5s2,s4 (23) 

and 

(c^,,_,,(ki)c^2,-s2(k2)cl,3,_,3(k3)cj;4^_,^(k4)) = 

(5(ki - k4)(5(k2 - 'k.z)^vx,vi^v.2,v3&sx,si^s2,si - 5(ki - k3)J(k2 - V^S^^^^^S^^^v^Ss^^s^bs^^si- (24) 

("2") 

As a result, contains four terms, associated with direct and exchange contributions to the signal at the Schrodinger 
detectors. 

In case both optical and matter- wave detectors are broadband [^, with response spectra centered at in and 
ujs respectively, Eq. (pO|) reduces to 



^(0 « '7m i '7s(l"3)%(r4) / ^7-4/ dT3 / dT2 \ dTi 

I Jo Ja Jo Ja 

X (E(-)(ri,Ti)E(-)(r2,r2)*t(r3,r3)*t(r4,r4)*(r4,T4)*(r3,r3)E(+)(r2,T2))E(+)(ri,ri)) 

rti rts rt2 rti 

+ 77s(r3,r4)rys(r4,r3) / dri / dr^ / dT2 / dri 
Jo Jo Jo Ja 

X (E(-)(ri,Ti)E(-)(r2,r2)*t(r3,r3)*t(r4,r4)*(r4,r3)*(r3,T4)EW(r2,T2))E(+)(ri,Ti)) 

nts nt2 Hi 

-77^(r3,r4) / (ir3 / dT2 I dri 
Jo Jo Jq 

X (E(-)(ri,ri)E(-)(r2,r2)*t(r3,r3)*t(r4,r3)*(r4,r3)*(r3,r3)E(+)(r2,T2)EW(ri,Ti))} (25) 
where the efficiency of the photodetectors 

77„(r,) = lim / dT{d{0)d{T))e'^"^^ (26) 

is assumed to be the same for all detectors. 

The "self-efficiency" rjs of the Schrodinger detectors is 



77.(r„r,) = hm ^ l],(r,)r!:(r,) f dr(vl>,(r,, 0)*]: (r„ r))e^(---=)- 



(27) 



with 

rjsir^r,) = Tjs{ri), (28) 

and their "cross-efficiency" is 

77^(r,,r,-) = 2 lim f drf dr' [ dr"e'('^^-"'=)(^+^') V e-'['^''i-i^+'^''2''2^'l 

'^"^■^0 Jo Jo ,^,,t^^,q, 

X {r.)\^\n,2 (r,)0q2 (r,)'/>q2 (r^) 

+ e-(----)(---")l):^(r,)f7,,(r,)l]:^(r,)17,,(rO|0q,(rOn0q,(r,)pl . (29) 
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In expression (]25|), the contribution of the Maxwell detectors is the same as in conventional photodetection theory, 
as expected since we use the same photoabsorption model. In contrast, the Schrodinger detectors are responsible 
for three distinct terms, resulting from the fermionic commutation relations of the electrons and ions. The first one, 
proportional to rig(r^)r]s{r4), is of the same nature as the expressions familiar from photodetection theory. It would 
take the same form in case the ionizing fields at and Vj lead to completely distinguishable signals. 

The second and third terms are of different nature, resulting from the quantum statistics of the Schrodinger 
detectors. The first of these, proportional to rig{vi, Vj)rjg(Yj^ Yi), originates in the exchange of the composite Schrodinger 
field ^i,(r), which is bosonic. We note that the dependence of this term on the product rjs{vi,Tj)rjs{Yj,Yi) rather 
than on an expression proportional to a fourth-order correlation function of ^^{v), is a direct consequence of the 
assumption that this field is initially in a vacuum. Finally, the third term, proportional to the cross-efficiency rj^^ 
involves the exchange of the electrons or of the ions only. Physically, it accounts for the fact that the electrons don't 
know their origin. This is the one term in Eq. (|2|) that depends explicitly on the fact that the composite bosons are 
annihilated and a pair of fermions is generated. The factor of 2 in rjx results from the fact that either the electrons or 
the ions can be exchanged. Hence, the first two terms can be thought of as being bosonic in nature, while the latter 
term is fermionic. 

We have previously mentioned that the (A^ 4- Af)-th order transition probability amplitude containsfiV -I- M)\ 
terms. Therefore 'P{t) is made of (4!)^ contributions, whose general form resembles that of expression (H). Hence, 
it is readily possible in principle to make measurements sensitive to both the bosonic and fermionic contributions to 
Vit). However, if one is only interested in the electron counting rate 

w(riii, . ..Viti) = —V{viti, . ..Vit4) (30) 

014 . . . Oti 

then the specific ordering of ii, . . . , i4 selects just one of these contributions. In addition, since w contains 4 partial 



derivatives, the fermionic contribution to Eq. (25) gives no contribution to the counting rate. Specifically, for the 
time ordering ii < t2 < ^3 < ^4 one finds readily 

w (ri^i, . ..Vit4) = ?7m(ri)77„j(r2) 

X {r/,(r3)ry,(r4)(E(-)(ri,ti)E(-)(r2,t2)*^(r3,t3)*^(r4,t4)*(r4,i4)*(r3,i3)E(+)(r2,t2))E(+)(ri,ti)) 

+ ?7,(r3, r4)77s(r4, r3)(E(-)(ri, ii)E(-)(r2, i2)*^(r3, i3)*^(r4, i3)*(r3, t4)EW (r2, t2)^^+\viM))} ■ (31) 

The first term in this expression is analogous to the standard fourth-order correlation function of optical coherence 
theory, with two of the electric fields replaced by the Schrodinger field. The second term, which is as we have 
seen an exchange term resulting from the undistingushability of the detectors, can be eliminated if they are rendered 
distinguishable. This could be realized experimentally in several ways, for instance by using ionizaing lasers of different 



frequencies. |20| In that case, the counting rate reduces to a simple measure of the normally ordered correlation function 



G(^+*'^)(xl,...,x^.+M) = 
where Xi = {vi,ti}. 

We conclude this section by noting once more that different gating schemes easily produce different cross-correlations 
functions between the optical and the matter-wave fields. It is easily shown that the gating scheme always selects the 
operators arrangement in such a way that time increases from the outside to the inside in the resulting correlations 
functions. As a result, of selecting the appropriate temporal ordering, it is possible to measure correlations functions 
involving different combinations of light and matter field operators. For example, the choice ^3 < ti < ^4 < t2 yields 
the correlation function 

(*t(r3,t3)E(-)(ri,ti)*t(r4,t4)E(-)(r2,t2)E(+)(r2,i2)*(r4,<4)E(+)(ri,ii)*(r3,i3)). (33) 



IV. EXAMPLE 



We now illustrate this theory by calculating selected third-order correlation functions of a system consisting of 
linearly coupled optical and matter- wave fields described by the effective Hamiltonian 
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JX o 

Ho = ^ fuvao-iaa + ^ hujiclci + (sf^alci + aPaaCi + /i-c.) , (34) 



where i? and S are the number of modes in the optical and matter field, respectively, with eigenfrequencies and 
LOi. In the following, we label the optical field with greek indices and the matter waves with roman letters. The mode 
annihilation operators 

a„ = yrf3ru;(r)E+(r) (35) 

and 

c, = y"dV0*(r)*(r) (36) 

are the bosonic optical and matter mode operators associated with the mode functions UQ,(r) and p/if j (r) , respectively, 
and gf'^ and g^f^ are coupling constants which allow for a parametric amplification type of matter-light coupling as 
well. 

The effective linear interaction of Eq. (^^ while not conserving the number of particles, is a good approximation to 
the description of the coupling (for instance via a two-photon Raman transition) between a macroscopically populated 
Bose-Einstein condensate ground state and some weakly populated side-modes, see e.g. [^,|3|. As long as the 
condensate is not significantly depleted by that interaction, one can invoke an "undepleted" mode approximation to 
replace the expectation value of the associated field mode operator by a c-number order parameter, leading to an 
interaction of the general form (|3j), with g^^'' and g^f^ proportional to the order parameter. 

In addition to their Hamiltonian evolution, we assume that the matter and light modes are also coupled to Markovian 
thermal reservoirs which result in exponential losses at rates and k^. The Langevin equations describing the 
evolution of this system are Q 

^x(<) = -Mx(t) + B^(i) (37) 
where ^(i) are reservoir operators with (^(t)} — and 

(38) 

In the general expression dH) x(i) = {aa{t),Ci{t)Y , a = {1, i = {1, S"}, and M is a (i? + S*) x (i? + 5) 

matrix whose elements depend on the system parameters, and the difusion matrix is D = BB^ = diag{naKa,niKi), 
where fia and are the thermal populations of the reservoir modes. 
We consider specifically the two-mode problem with drift matrix 

M=(^]/^ (39) 



D = I 1 _ " ) (40) 



and diffusion matrix 

niK2 

and evaluate the stationary value of the third-order correlation functions 

Gi^) = (*t(i.,,t,)*t(i.2,t2)E(-)(r3,i3)E(+)(r3,t3)*(r2,i2)*(ri,ii)) (41) 

and 

G(3) = (*t(i.^^t^)E(-)(r3,t2)*t(r2,t3)*(r2,i3)E(+)(r3,t2)*(ri,ti)), (42) 
using the derivation outlined in Appendix A. 

(3^ 

The normalized stationary cross-correlation function gx {ti,T2) is shown in Fig. 1 as a function of the time 
differences ti = — ti and T2 = ^2 — ^i- It corresponds to the time gating ti < t2 < t^, so that T2 < ti. For the 
particular parameters at hand, and for a fixed time difference ti we observe a bunching-like behavior as a function 
of the time delay T2 between the Schrodinger detectors. This should be contrasted to Fig. 2, which shows the 
normalized correlation function gy (ri, T2) for the same parameters as in Fig. 1. In this case, T2 is the time difference 
between the photodetector and one of the Schrodinger detectors, hence the behavior of the correlation function along 
T2 exhibits a different behavior for fixed ti. Note that the line ri — T2 yields identical information in both cases, since 
it corresponds to the identical physical situation where the optical detector and one of the matter-wave detectors are 
simulataneously turned off. 
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V. SUMMARY AND CONCLUSIONS 



The dynamical interplay between Maxwell and Schrodinger waves, which is the cornerstone of most novel appli- 
cations of nonlinear and quantum atom optics, requires one to be able to manipulate and characterize not just the 
statistical properties of the optical and matter fields individually, but also their cross-coherence properties. 

In this paper, we have extended the theory of mutual coherence between optical and matter waves to the study 
of higher-order correlation functions. While in most studies of degenerate atomic quantum gases it is sufhcient to 
treat them as (composite) boson or fermion fields, this is no longer the case if the detection mechanism destroys the 
particles, as is the case for photoionization. The quantum statistics of the constituent particles play an important 
role, as illustrated by Eqs. ( ^7|) and ( p9[ ) which show a bosonic term involving the simultaneous exchange of ions 
and electrons, and a fermionic term involving the exchange of ions or electrons alone. While it is possible to conceive 
specific counting schemes insensitive to these terms, it would be extremely interesting to observe them. 

The generalization of this analysis to fermionic fields is also of much interest, since in this case, the constituents 
involved are a fermionic field and a bosonic field instead of two fermionic fields as was the case in this paper. The 
exchange contributions to the ionization probability will be quite different in that case, since electron exchange, as 
well as the simultaneous exchange of ions and electrons, wilU now lead to fermionic contributions, while the exchange 
of ions produces a bosonic term. 
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APPENDIX A: CORRELATION FUNCTIONS MATRIX ELEMENTS 

From the expansions (|3^) and (^6|) of the Schrodinger field operators, the correlation function G'i'^^ is 
Gl'^ = E 0*(ri)0*(r2)<(r3)u^(r3)0fc(r2)0Kri) 

ijal3kl 

{ (f/-^a(t3,ti)(t/-^)t)^,(c/-^e(t2,t2)(c/-^)0,^.(c/-^a(ti,i3)(^-^)t),„ 

+ (f/-^e(t3,tl)(t/-^)t)^,(C/-^e(tl,t2)(C/-^)0y(^"'^(^2,i3)(^-^)0fc« 

+ (f/-^a(t2,ti)(t/-^)t),^(c/-ie(t3,t2)(c/-^)t)^^.(c/-ia(ti,t3)(^-^)t),„ 
+ {u-'g{t2,h){u-y)^^{u-'g{tuh){u-')^)^^{u-'g{t,,t,){u-y)^^ 
+ {u-'g{h,h){u-y)^^ {u-'g{ts,t2){u~y)^^ {u~'g{t,,ts){u~y)^^ 
+ {u-'g{h,h){u-y)^,^ {u-'g{t,,t2){u-y)^^ {u-'g{t3,ts){u-'y)^j (ai) 

where U is the matrix of eigenvectors of the system matrix M, so that 

U-'^MU = diag{Xk), k = {1, R + 5}, (A2) 
where {Xk} are the eigenvalues of M, and 

g,, {t, s) = {UBU% ^^^^ t > s, 



g,^{t,s) = iUT)U%^^^^ t<s. (A3) 



'3 



(3) 

The correlation Gy has a similar expression. 
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The six terms in the Eq. ( |Al[ ) follow directly from the steady state solution of Eq. (^7|) given by 

x(t) = f e-^'^'-*">B^{t')dt'. (A4) 

J —oo 
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FIG. 1. Normalized correlation function {ti, T2) as a function of n = ts — t\ and T2 = t2 — ti for ki = 0.15, K2 = 0.25, 
ni = 0.01 and n2 = 0.1, in units of the coupling constant g. 
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